Ferromagnetic resonance with a magnetic Josephson junction 
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We show experimentally and theoretically that there is a coupling via the Aharonov-Bohm phase 
between the order parameter of a ferromagnet and a singlet, s-wave, Josephson supercurrent. We 
have investigated the possibility of measuring the dispersion of such spin waves by varying the 
magnetic field applied in the plane of the junction and demonstrated the electromagnetic nature 
of the coupling by the observation of magnetic resonance side-bands to microwave induced Shapiro 
steps. 



Rotation symmetry associated with the 0(3) orthog- 
onal group forbids a coupling between a scalar s-wave 
order parameter and the vector order parameter M of 
the ferromagnet. However a Josephson junction [l| de- 
fines a plane, the 0(3) symmetry is broken, and such a 
coupling is possible. Here we describe a part of the rich 
spectroscopic magnetic resonance possibilities that this 
observation implies. It is possible to perform a "photon 
free" FMR experiment [|| on about 10 7 Ni atoms, some- 
thing infeasible with standard FMR techniques. 

Interactions in nature reflect certain gauge groups and 
the associate phases which generate a vector potential 
A, called the Berry connection^. Interactions via elec- 
tromagnetic fields generated by the conserved electrical 
currents reflect the U(l) gauge group and the Aharonov- 
Bohm (AB) phase[H. Associated with angular moment 
is the gauge symmetry SU(2) and the familiar Lie al- 
gebra, i.e., the spin commutation rules. The AB is re- 
placed by the spin Berry phase. It is often imagined 
that magnetic moments might interact with the Joseph- 
son current by direct spin- flips, which would involve the 
spin Berry phaseJBJ. In this paper it will be shown 
that quantitatively the interaction between the Joseph- 
son current and the order parameter in superconduc- 
tor/fcrromagnet /superconductor (SFS) junctions can be 
explained in terms of the AB-phase and regular electro- 
dynamics. It will be shown such an experiment measures 
rather directly the magnet correlation function. 

We fabricated Nb/Pdo.gNi .i/Nb Josephson junctions 
by in-situ angle evaporation through a resist mask 
and subsequent lift-off. The mask is defined by e- 
beam lithography on a Polyether Sulphone - PES (500 
nm)/Si 3 N 4 (60nm)/Polymcthyl Methacrylate - PMMA 
(350 nm) trilayer [fj , and etched in a Reactive Ion Etch- 
ing (RIE) chamber. The SiaN4 is etched 1 minute 30 
seconds in SF§ and the PES in Oxygen plasma for 10 
minutes, giving an undercut of 500 nm. The mask fabri- 
cation process is schematically presented in Fig. QJ,. The 
SisN4 suspended bridge allows shadow evaporation. A 
scanning- electron-microscope (SEM) picture of the mask 



including the suspended bridge is reported in Fig. [TJi. 
The first Nb layer is evaporated at -45 degrees with re- 
spect to z axis while the PdNi is evaporated at 45 degrees 
and the second Nb layer at 47 degrees. The overlap in the 
y direction defines the junction area. The shadow evapo- 
ration is illustrated in the lower drawing array of Fig. [TJi. 
Eight junctions are evaporated on the sample chip. A 
SEM image of one of the junctions after lift-off is shown 
in Fig. [T]i. The electron-gun evaporation is carried out 
in ultra-high- vacuum (UHV) with a base pressure lower 
than 10 -9 mbar. The two Nb superconducting films are 
50 nm thick and the ferromagnetic layer of Pdo.gNio.i 
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FIG. 1: Description of the sample, (a) The fabrication proce- 
dure. Trilayer mask is etched to obtain the suspended bridge 
and the undercut, then the junction is fabricated by shadow 
evaporation and liftoff, (b) The principle of the experiment. 
AC Josephson current across the junction excites the spin 
wave modes which in turn couple to the Josephson current and 
rectification takes place, (c) Scanning Electron Microscopy 
(SEM) picture of the SiaN4/PES mask before angle evapora- 
tion, (d) SEM picture of a junction after lift off. 
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is 20 nm thick. The Ni concentration is measured by 
Rutherford Backscattering (RBS) on a test sample. The 
magnetization loops obtained by SQUID magnetometry 
with in-plane and out-of-plane magnetic field show a pre- 
dominant perpendicular anisotropy. Finally, a schematic 
view of the junction including the principle of the exper- 
iment is to be found in Fig. [T|d. 

Typical current- voltage (/V) characteristics are shown 
in Fig. [2] as function of temperature. Well below the 
critical temperature, the critical current I c w 7 /iA and 
normal resistance R n = 1 give a Joscphson coupling 
of 7 /iV, consistent with early studies on highly under- 
damped PdNi-based Josephson junctions The IV 
characteristics are not hysteretic confirming overdamped 
phase dynamics and are well described by the resistively- 
shunted-junction (RSJ) model pj. The critical current 
versus temperature (see inset of Fig. [2]) shows the typ- 
ical linear behavior expected when the Thouless energy 
of the ferromagnetic layer is larger than the Nb super- 
conducting energy gap. This linear dependence has been 
observed previously in highly underdamped j unctions Q. 
The junction critical temperature is about 7.0 K while 
the critical temperature of the Nb leads is 7.6 K and 
their critical current over 500 /^A at low temperature. We 
have measured four ferromagnetic junctions on the same 
wafer, the dispersion of the critical current from junc- 
tion to junction is about ±1 /iA, AR n is 0.15 while the 
I c R n varies by less than 3% from junction to junction. 
We have also fabricated non-magnetic junctions by the 
same process but replacing the PdNi with a thicker 70 
nm Pd layer. These junctions have a much larger critical 
current of about 44 /iA. 
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FIG. 2: Current-voltage characteristics with increasing tem- 
perature. The curves have been shifted vertically for clarity. 
They do not show any hysteresis as expected when the phase 
dynamics is strongly damped. The critical current as a func- 
tion of the temperature is shown in the insert, the dashed line 
is a linear fit. 



For a square junction of side L, the total super-current 
is given by the integral [l[ 

,L/2 ,L/2 

I S = J C / dx dy sin <f>(x, y, t) (1) 

J-L/2 J-L/2 

with 

0(r, t) = fa + ujt - ^ J A-dr, (2) 

where <po is an arbitrary phase, wj = (2e/h)Vo is the 
Josephson frequency, and the last term is the AB phase 
Q, involving the vector potential A. We use a gauge 
where A = A{f, t) z, the direction i being perpendicular 
to the junction surface. Therefore <p(f, t) = 4>o + kx + 
ujjt — 4>i, where <f>i reflects time dependent fields and 
k = [Aed/K)iiQH + (4ea/h)fj,oMo y . Here Mq v is the y- 
componcnt of the static magnetisation Mo, the applied 
field H is in the y direction, 2a and 2d— 2 (a + A) are the 
actual and magnetic thickness of barrier and A the Lon- 
don penetration depth. These equations describe both 
the statics and the dynamics of our junctions. 

The approach is similar to that used for junctions with 
magnetic impurities in a normal metal barrier It is 
necessary to determine the appropriate boundary condi- 
tions for the solutions of Maxwell's equations. For rea- 
sons of transparency, it is not at all useful to solve the 
very difficult problem in which the solution within the 
barrier is matched to that in the exterior region to the 
junction. Within the junction we can ignore the dis- 
placement and transport currents since the wavelength 
of light A and the skin depth 6 are both larger than the 
dimensions of the junction at the Josephson frequency ojj 
relevant for the FMR. We observe that the impedance of 
the junction of ~ 1 fl is much smaller than 377 f2 of free 
space and as a consequence there is essentially no radi- 
ation from the junction. The displacement current, and 
evidently, the transport current can therefore also be ig- 
nored in the exterior region. It is therefore only necessary 
to integrate Ampere's circulation law 

Vx#=J 8 (r, t)i- J s = J c sin cf>(r, t) (3) 

where <f>o is absorbed by a time translation, and <f>(r, t) = 
kx + LJjt — 0i . Required is the additional AB phase shift 

with the magnetic system reflected in A 1 . In linear re- 
sponse 0i is considered as a perturbation and the dc sig- 
nal is 

I1= -piu T dt r dx r dy ^ Al (5) 

UwjTJq J_ l/2 J_ l/2 dt 
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which includes a time average over a single period T. The 
determination of A 1 requires first the vector integration 
of V x H = J s (r,t)i and then V x A = B with B = 

Ho(H + M). Even with the simplifications of the previous 
paragraph, this is an involved calculation. It is useful to 
make some formal manipulations in order to avoid this 
double integration. First Ii is written as 

4ae i r T r L/2 r L/2 f o ~ -\ - 

(6) 

using Ampere's law J = V x H . Performing an integra- 
tion by parts on time we have 

4ae 1 f T , f L/2 , f L/2 , /-. f+\ OA 

m 

Using the fact that the Poynting vector, and hence H x 
dA/dt — on the surface, this is integrated again by 
parts using V • (H X (dA/dt)) = (dA/dt) • V x H -H ■ 
V x (dA/dt), to give 

h = ^ / dx dyH-^-yxA. (8) 

nivjT J J_ L/2 J_ L/2 dt 

Then, since Vxi = B and B = ^o(H + M), the signal 



where dv = 2adxdy is the elementary volume, the inte- 
gral is over the volume of the magnetic layer, and the 
average is indicated by the bar. The resonance of the 
ferromagnetic layer is contained in x 2 (i), the dynamic 
susceptibility, and 



M t (f,t) = Id? I dt l Xl (r-r;t-t , )H l (^,t') 1 (10) 

where i = x,y,z. The simple expressions Eqs. ((9]) and 
(fTU|) are a principal theoretical result presented here and 
have an obvious interpretation in terms of the magnetic 
energy They demonstrate that Joscphson junction mag- 
netic spectroscopy measures very directly the magnetic 
susceptibility correlation function Xi(r — r*;t — t'), and 
all the other excitations to which that couples, in much 
the same manner as does neutron scattering. As will be 
seen below, the advantage is that this technique couples 
preferentially to small q wave vectors. 

Since Hi(r, t') is periodic in time, the time convolution 
Eq. (JTUJ), reduces to a product and if the susceptibility 
is sufficiently local and then 



Mi(r,wj) = (x^cjj) + itf(uj))Hi(?,wj) (11) 

in the usual complex notion. If rather the response is 
non- local then 

Mi(k,uj) - (XiiKu^ + iXi&ujMiiKujj) (12) 

and it requires a Fourier expansion of the spatially de- 
pendent Hi(r,u)j). Finally, when Uj ~ I C R, as for the 
lowest voltage experimental signals, the linear response 
approximation is not strictly applicable and high har- 
monics of loj must be accounted for. Similar expressions 
apply but now, in particular, "half-harmonic" signals oc- 
cur since the super-current contains higher harmonics 
and can, corresponding to the second harmonic, excite 
a resonance at uj s when lo,j/2 = u) s . 

Now all that is required is a single integration of 

V x H = J c sinffcx + ujt]z, (13) 

but which is not a simple task. It is trivial to ver- 
ify by differentiation that such an integral is H p = 
— ( J c /k) cos(kx + ujjt)y. However this does not satisfy 
the boundary condition that the current density is zero 
outside the square junction region. With the present 
gauge A = A z z, it is necessary to find a solution of 
Laplace's equation V 2 A Z = such that H = H p + Hi 
where Hi = V x A is such that H satisfies Eq. (|13[) inside 
the square but has V x H = outside. A little reflection 
suggests there are two contributions to Hi which must 
be accounted for. First, in general, reflecting the even 
part J c cos kx sin ujjt of the current density there is a net 
oscillating super-current I C (H) which causes a circulat- 
ing magnetic field about z, and second, associated with 
the spatially odd part of J c s'm[kx + uijt], there is a uni- 
form component of the field in the y-direction. The first 
contribution is determined by considering the problem 
with k = 0, i.e., with a uniform current density J s . The 
solution is H s = ^J s (xy — yx) cos w jt + . . . where the el- 
lipsis reflects the relatively small corrections for a square 
as compared with a circular cross section. In what fol- 
lows this correction is ignored. The corresponding vector 
potential has A z = (1/4) J s (x 2 + y 2 ). For finite H, inte- 
grating the even part of the current density gives an aver- 
age super-current density of J s = J c sin(kL/2)/(kL/2). 
By inspection it is observed A z = (l/A)J s (—x 2 + y 2 ) 
satisfies V 2 A Z = 0. The corresponding odd H° = 
xy — yx) cosujjt + The sum H p + H° correctly re- 
duces to H s in the limit k — > 0. It is the case that H p 
(and H p + Hf) implies a uniform oscillating field but one 
which diverges as k — > 0, whereas physically, the even 
part of H p should be proportional to k reflecting the 
Josephson screening of fields. That the tangential ap- 
plied field H be continuous requires the current, induced 
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by the even part of H, to be zero at the surface. The 
even part of H p is —J c cos k kx sinujjty(coskx/k) sinujty 
The required even part of Hi is now Hf = 

J c cos , fc 2 sina;jty(cos k^/k) smujjty, which is equally di- 
vergent as k — > 0. The net result of integrating Eq. (fT3|) 
is therefore 

-» J kL 
H = -^[cos(fca; + ujjt) — cos — sinwjtjy 



+ o Jc fcL (-^y-2/x)cOSCJ,/^ + ... (14) 

Imagine that the magnetic layer is composed of a num- 
ber of independent crystallites so that the response is 
local and Eq. (fTT]) would apply. This local assumption 
also has the merit of being an useful illustration of the 
theory since it leads to a relatively simple prediction of 
the H dependence of the signal which can be compared 
with experiment. This helps determine if the response is 
indeed local, or extended. There are some complicated 
integrals involved in the evaluation of Eq. (|9|) . The result 
is written, in closed form, as 




where $ r f = (2aL)B r f = (2aL)/j Ic(0)/L is the flux due 
to the radio frequency field and where 



and 
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FIG. 3: In blue is F y and in pink F x . These curves satisfy 
the evident requirement that the coefficient of the x and y 
responses be the same when B = 0. Markers denote the 
experimental values of the resonance amplitude as function 
of the applied field, with the normalization constant as the 
fitting parameter. 



with x = kL, reflect the geometrical structure of the 
coupling. The equilibrium magnetization is along the z 
axis, and the magnetic resonance signal is contained in 
x"(a>j) and x^'(wj), the Fourier transforms of the imag- 
inary part of the susceptibility. The two functions F x 
and F y are plotted in Fig. [3] As required by symmetry 
F x = F y when k = 0. The x"(wj) response near the first 
zero in I c (Bo) = I c is about four times the maximum re- 
sponse to Xy(u>j). While there is a clear reflection of the 
Fraunhofcr diffraction pattern in the Xy{uJj) response, 
the l/k 2 , i.e., 1/H 2 response dominates that to x"(wj) 
with only modulations due to diffraction effects. With a 
single flux quantum threading thorough the junction the 
F x is zero reflecting the net absence of a circulating cur- 
rent. In contrast F y is a maximum since the current is 
odd and the junction constitutes a small flat solenoid car- 
rying the critical current density and hence a maximum 
field internal to the junction. 

Given a static magnetisation M — M z z the magnetic 
susceptibility might be approximated by 

X»(fc, u,) = xgft a,) » 7e M, g K + % -p u)i (16) 

where r is the relaxation time, u> s the frequency of the 
FMR mode and a the spin-wave stiffness. Here j e = 
jUg/Ti with /is the Bohr magneton. 

In Fig. 2] we report the dynamical conductance at zero 
applied magnetic field (solid line) and the theoretical 
expectation (dotted line) from Eq. (p~5|) and (fTBf with 
a = 0, using uj s and 1/t as a parameters. We ob- 
tained u s = 23 /iV as expected fro m the Kittel's formula 
uj s = "/ c \/{Hk — ^ttMs) 2 — H 2 , where the anisotropy 
field Hk and saturation magnetisation Ms have been 
measured separately by SQUID magnetometry The 
FMR frequency is consistent with the value obtained in 
a reference sample by EPR spectroscopy Q, i.e., the ob- 
served signals are entirely consistent with the coupling of 
the FMR resonance to the superconductivity via the AB- 
phase. The resonance at 11.5 /xV is a subharmonic of the 
main mode. Kittel's formula predicts that uj s decreases 
with increasing H whereas the half-harmonic signal of 
Fig. @] actually seems to increase (see inset of Fig. 2] ). 
This leads us to believe that this signal corresponds to a 
finite value of a and k in Eq. (|16[) . That the spin- wave 
dispersion might be measured in our type of experiment 
is an exciting possibility warranting further investigation. 
Moreover, in Fig. [3] we also report the amplitude of the 
FMR signal (markers) as a function of the applied mag- 
netic field. The experimental data follow the F x coupling 
function obtained above. 

The electromagnetic nature of the response can be con- 
firmed by a study of Shapiro steps [H, Gl|. In the absence 
of a magnetic resonance mode, an applied radio frequency 
field gives rise to such steps. In order to account for these 
we write for the bias voltage across the junction 

V = Vo + v cos fit 
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where v and SI are the amplitude and frequency of the 
applied microwave field while the constant Vq, as usual, 
corresponds to a Joscphson frequency uij = 2eVo/Ti. The 
Josephson current is now J c times 

2ev 

sin[fca; + cut + — — sin Sit]. 



Expanding this sine with the Jacobi- Anger identity gives 




with S2n = X)± sin(fcz + ujot + 2nSlt) and S2n-i = 
± sin(A;a; ± ojot + (2n — l)Slt), where So = S , 2n=o- 
Each term is of the form 



J™ sin[fcz + u> t ± nQi] -> J™ sin[fcz + (17) 

where, in well known fashion[l|, J™ = J c J n (2ev/1iSi) in- 
volves the appropriate Bessel function J n (2ev/TiSl). For 
the voltage at which u n = there is a direct contri- 
bution J™ sin (0o + kx) to the average current which is 
equivalent to the zero voltage critical current step but 
displaced to V n = nhSl/2e and reduced from I C (H) to 
I c (H)J n (2ev/hSi). This corresponds to the principal 
Shapiro steps, shown on the top panel of Fig. [5] For 
the magnetic response, in the linear response regime, the 
theory developed without an applied microwave field can 
be adopted. All that is needed is to replace J c by J™ and 
ujj with cj" in the appropriate places as described above. 

We have investigated the Shapiro steps for different 
microwave power. The top part of Fig. [S] shows the ap- 
pearance of the Shapiro steps in the current-bias char- 
acteristics for increasing microwave power. The ampli- 
tude of each step follows the appropriate Bessel func- 
tion J n (2ev/hSl), as expected. In the bottom panel of 




FIG. 5: Top panel: IV curves showing the Shapiro steps, 
taken at T=1.3 K and with the microwave frequency of 1 GHz 
with different incident power. Bottom panel: the dynamical 
resistance showing Shapiro steps and side-band resonances, 
taken at T=35 mK and with microwave frequency of 17.35 
GHz. In passing from the bottom to top the power increases 
from dBm to 20 dBm, in steps of 5 dBm. 



2.0 



t51.5 



1.0 




10 20 30 40 50 

V UlV) 



FIG. 4: Dynamical resistance of a ferromagnetic Josephson 
junction showing ferromagnetic resonances with first and sec- 
ond harmonic. Insert: the magnetic field dependence of the 
resonance at the second harmonic. The fit takes into account 
the spatial dispersion of the spin wave mode, see text. 



Fig. [5] is presented the dependence of the magnetiza- 
tion induced side-bands on the micro- wave power. To 
make these bands more evident, shown is the dynami- 
cal resistance as a function of the voltage. It is clear 
from the data that the amplitude of the side-band reso- 
nances follows the Bessel function of the Shapiro steps, 
as expected from the theory described above. Indeed 
the Shapiro steps can be seen as "replica" of the critical 
current at finite bias and hence the side-band amplitude 
follows that of the steps. The two sidebands correspond 
to the two poles of the dynamical susceptibility, Eq. (|16[) . 

In conclusion, we have described experiments and de- 
veloped theory to demonstrate that the relative AB phase 
of the superconductors which comprise a Josephson junc- 
tion couples to the magnetic order parameter of a ferro- 
magnet. We have thereby performed an FMR experiment 
with a sensitivity which greatly exceeds that of conven- 
tional cavity FMR. Since the coupling is via the magnetic 
field it is not necessary to have the current pass through 
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the magnetic material. It might be imagined that the 
magnetic layer be the top layer of a FSIS structure in 
which the adjacent S-layer has a thickness of the order 
of, or less than, the London penetration length. The 
possibility of measuring the dispersion of spin- wave exci- 



tations has also be investigated. Our method^, @> El of 
coupling superconductivity to magnetism measures di- 
rectly the dynamic susceptibility x"(<f, w) with an en- 
hanced sensitivity for small wave- vectors, complementary 
to neutron scattering. 
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